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New data on neutrino masses is very exciting. It is the first evidence for physics 
beyond the standard model. In this talk we discuss this data and how it may 
(or may not) be used to constrain theories of charged fermion masses. The talk 
is divided into two parts. In part I, we briefly discuss the theoretical framework 
for fermion masses and mixing angles as well as some of the most salient data for 
neutrino oscillations. In part II, we define the framework of a "predictive" theory 
of fermion masses. We then consider a particular theory which makes significant 
testable predictions. 

1 Introduction to Neutrino Oscillations 

1.1 Theory 

Consider the fermion mass sector of the minimal supersymmetric standard 
model[MSSM]. Fermion masses are defined in the superpotential \W]. Three 
3x3 complex Yukawa matrices A^, Ad, Ag are needed to describe up, down and 
charged lepton masses, while is needed for neutrino masses. In the effective 
low energy theory, the first three matrices contain 13 observable parameters 
- 9 masses and 4 quark mixing angles. The superspace potential is given by 



where all fields are described by two-component left-handed Weyl spinors. 
The two Higgs doublets iJ„, H^l couple to the quark and lepton doublets 



and the left-handed anti-quark u, d and anti-lepton e, v singlets. Finally, 
we have added three sterile neutrinos u into the MSSM in order to describe 
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neutrino masses and mixing angles. They are sterile since they have no elec- 
troweak quantum numbers. As a result, a 3 x 3 Majorana mass term with 
mass matrix M may be added to the theory without breaking any local gauge 
symmetry. When the eigenvalues of M are taken much larger than the weak 
scale we obtain the well-known see-saw mechanism.l!l We return to neutrino 
masses shortly, but first consider the charged fermion masses. 

In a fundamental theory, presumably defined at some large scale, the 
quark and lepton Yukawa matrices are complex 3x3 matrices which encode 
fermion masses and mixing angles. In this basis, labeled by the subscript 
0, the Yukawa matrices are not diagonal; while electroweak interactions are 
diagonal. This is the so-called weak eigenstate basis. At the weak scale 
the electrically neutral Higgs components obtain vacuum expectation 

values[vev] ^ = ^ = respectively, where v = {^/2Gf)-^/^ = 

246 GeV and Gf is the Fermi constant. We then identify the charged fermion 
mass matrices rriu = \u and similarly for down and charged lepton masses. 
The u, d, e mass matrices can each be diagonalized by two specific unitary 
matrices. For example, for up quarks we have 



rric 

mt 



= Ul.^muUiiu- (3) 



In this mass eigenstate basis, electroweak interactions of quarks are no longer 
flavor diagonal. Moreover, the CKM matrix is given by the mismatch between 
diagonalizing left-handed up and down quarks. We have 

VcKM = C/l„ ULd (4) 

Now consider leptons. The charged lepton mass is given by 

me = Ae^ (5) 

and is diagonalized by the transformation Zq = ^ C^le and similarly for left- 
handed anti-leptons such that 

eomeeo = e me C/_Re e (6) 




For neutrinos we now have 

lUl^KDoHu MDo (7) 
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Figure 1. Diagram generating the effective neutrino mass operator 



The mass scale M may be much larger than the weak scale; thus in the limit 
M » Mz we can integrate the sterile neutrinos Do out of the superspace 
potential and obtain the effective dimension five operator, see figure 1 - 



{lUl^KH^)M-HHuXlUlJ) (8) 

This is the Gell-Mann, Ramond, Slansky and YanagidaEI seesaw mechanism 
for obtaining light Majorana neutrinos. 

The dimension-5 operator is defined at the scale M and must be rcnor- 
malized down to Mz where obtains a vev. The Dirac neutrino mass is 
defined by nii, = and finally the effective Majorana neutrino mass is 

given by 



m 



l^f =Ul,m,M-^mlUl, (9) 




It is defined in the basis in which charged lepton masses are diagonal. This is 
the so-called flavor basis for neutrinos. It is in this basis that a W boson takes 
an electron into an electron neutrino, etc. ml^^ however is not necessarily 
diagonal. It can be diagonalizcd by a unitary transformation such that 



= C/t miff U* (10) 



The unitary matrix Uai and mass eigenvalues are the observables in 
neutrino oscillation experiments. The indices a = {e, fi, t} label the neutrino 
flavors, while i — {1, 2, 3} label the mass eigenstates. 

So consider a process in which lepton a turns (via weak W exchange) 
into neutrino flavor Va, a linear superposition '^^Uiai'i(0) of neutrino mass 
eigenstates, at position 0. A distance L away the neutrino ^{{L) = 
exp{imf L /2E)ui{0) is now detected as a linear superposition of neutrinos 
of flavor (3, see figure 2. This is the process of vacuum oscillation and the 
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Figure 2. Neutrino of type o produced at point and detected as neutrino of type f3 a 
distance L away 



probability for the state i>a to remain Vq, is given by the formula 



P{Va -> I^q) = 1 - 4X)fc<, PajPak sin^ A^fe 



where 



A,fe =1.27 



eV^ / {E/MeV) 



(11) 



EE|C/„,f (12) 

' ^'^%\ {L/ meters) 



Note, when neutrinos travel through matter there are possible matter 
enhanced neutrino oscillations which were discovered by Mikheyev-Smirnov 
and Wolfenstein [MSW] .□ These effects would take more time to describe and 
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are beyond the scope of this talk. In what follows, we shall specifically state 
when MSW oscillations are the dominant effect. 



1.2 Neutrino Oscillation Data 

We briefly discuss the evidence for atmosphcriJi, solarS and liquid scintillation 
neutrino detector [LSNDJeI neutrino oscillations. 

Atmospheric neutrino oscillation data 

Cosmic rays incident on the upper atmosphere scatter on heavy nuclei 
and produce neutrinos. The neutrinos are predominantly by-products of pion 
and kaon decays which result in two muon type neutrinos for each electron 
neutrino. These neutrinos are detected at Super-Kamiokande [SuperKju, a 
large water cerenkov detector located in a deep undergound laboratory in 
Japan. Electron and muon neutrinos scatter in the water producing electrons 
and muons which are distinguished by their cerenkov rings. Muons produce a 
sharp well-defined ring, while electrons multiple scatter in the water producing 
a ring with fuzzy boundaries. These atmospheric neutrinos are produced 
uniformly around the globe and the distance they travel to the detector is 
a function of the zenith angle ijj. Neutrinos coming from above {ip = 0) 
typically travel a distance L ^ 10 km, while neutrinos from below the detector 
(-0 = tt) must travel through the earth corresponding to L ~ 10, 000 km. Thus 
atmospheric neutrinos are a natural source for testing neutrino oscillations. 

Experimentalists typically analyze their data in terms of a two neutrino 
mixing model with mixing matrix 

U - -( ^"""^ (13) 
\- sine cose ) ^^^> 

For muon neutrino oscillations a = {/^, x} where x can be either e, t or s (a 
sterile neutrino). The persistence probability for is given by 

^K^^^M)-l-sin'20sin2(^) (14) 

where the muon neutrino energy E is inferred by measuring the energy of the 
muon it produces. _ 

The SuperK data u show (see figure 3) 

• # depends on the zenith angle V' or equivalently on 

• # z^e is independent of L/E. 
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Figure 3. Super-Kamiokande data for electron and muon neutrinos divided by Monte Carlo 
generated data assuming no oscillations plotted vs. L/E(km/GeV). 



This simple observation suggests that muon neutrinos oscillate while electi' 
neutrinos do not; hence I'rr 7^ i^e- After 535 days of running the best fi 

v-r) oscillations (including the fully contained, 



see figure 4d, assuming (i^^, 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 6 



partially contained and upward going muon data) is given byQ 

sm^29 =1.05 (15) 
Sm^ = 2.2 X 10-3 

and x^/d-o.f. « 64.8/67. We'll come back to the question of whether or not it 
makes sense to consider unphysical values of sin^20 > 1. If instead (z^^ — > Ve) 
oscillations are assumed, the fit becomes worse with x^/d-O-f ~ 87.8/67. 

In fact, there is additional experimental data which also suggests that 
j/j. Ve- This comes from the Chooz experimentQ, a neutrino detector which 
sits about 1 km from a nuclear reactor in France. The reactor produces 
electron neutrinos with a known flux. The experiment looks for the disap- 
pearance of these neutrinos. It is sensitive to the region of parameter space 
with Sm^ > 10~'^eV^ and sin^29 ~ 1. Their null result suggests that muon 
neutrinos cannot oscillate into electron neutrinos with these parameters, since 
otherwise they would have observed a disappearance of electron neutrinos due 
to their oscillation into muon neutrinos. 

Thus it is safe to conclude that i>x ^ for atmospheric neutrino oscil- 
lations. Hence we have two alternatives, v,j, — v^^ or Vg. How can this be 
reconciled? 

We consider a recent analysis by K. Scholberg at SuperK.Ll At the moment 
both possibilities fit the zenith angle dependence well. It should however be 
possible with better statistics to distinguish vs. by the zenith angle 
dependence. This is because there is an MSW effect in the earth for Vs but 
not for Vr- 

Also by considering the ratio of neutral current [NC] to charged current 
[CC] processes one can distinguish between the two. In this case there is 
preliminary data which favors = Vr- This ratio satisfies 

R{NC/cc) < 1 for ly^ = Us (16) 
= 1 for = J^r- 

Using SuperK data for tt" events produced by neutral current neutrino scat- 
tering in the detector one measures (see figure 5) 



^iNC/CC) = (^u/e)Mltc.Wo (1^) 

= 1.11 ± 0.06(data stat.) ± 0.02(MC stat.) ± 0.26(sys.) 



"More recent SuperK datalZI after 736 days of fully contained events and 685 days of partially 
contained events now gives a best fit value of Sm'^ = 3.2 X IQ-^eV^and sin 29 = 1.05 
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Figure 4. Super-Kamiokande fit to atmospheric neutrino oscillations assuming i/p — » v-r- 
The figure is taken from SuperK 98. The large dot corresponds to the result after ^ 700 
days of running. 



Solar neutrino oscillation data 



The flux_of solar electron neutrinos from the sun can be calculated, see 
for exampleJj The result depends on the details of the solar model, but in- 
dependent of which solar model is employed one finds disagreement with the 
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Figure 5. Neutral current process resulting in a tt — > 27 where the subscript a refers to 
any active neutrino. 



Table 1. Solar neutrino data 





ResultH/Theorja 


Homestake 

Kamiokande 

GALLEX 

SAGE 

SupcrK 


.33 ± .029 
.54± .07 
.60 ± .06 
.52 ± .06 

.474 ± .020 



data (see table 0) .| 

Solar neutrino data is typically analyzed in terms of two neutrino oscilla- 
tions where Ve — > with a; = a, an active neutrino {v^, Vr) or x = s, a sterile 
neutrino. The results of many fits to all the data give three possible solutions 
with X — a and one solution with x = s (see table These solutions fall into 
three categories — (1) small mixing angle [SMA] MSW, (2) large mixing angle 
[LMA] MSW and (3) the so-called " Just so" vacuum oscillation solution. 

LSND oscillation data 



''Note the most uncertain component of the neutrino flux is the highest energy so-called 
hep neutrinos. This affects the prediction for the high energy tail of solar neutrinos. 
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Table 2. Solar neutrino fits 



MSW 


SMA sin^20 ~ 10"^ - 10"^ 
= Us, Ua &m? ~ 10-'5 - 10-5 eV^ 




LMA sm^26»~.6-.9 
= Va Sm? 10-5 - 10-4 eV2 


Vacuum 


"Just so" sm^26» ~ 1 

= Ua Sm? ~ 10"" - 10^10 eV^ 



The LSND experiments uses a proton beam at Los Alamos National Lab 
to produce pions in a beam dump. Both tt+s and ■k~s are produced. In their 
stopped pion data, tt+s decay at rest into i-i^i^fx and /i+ — > Ve v^^. Note, 
there are no VeS in this decay chain. ■n~s, on the otherhand will typically be 
absorbed into a nucleus before they have a chance to decay. 

They then look for the signature of Ve in the detector, see figure 6. The 
Df, scatters on a proton and produces e+ n in the detector. The positron is 
observed via its cerenkov radiation. A short time later the neutron finds a 
nucleus and is captured, producing a gamma ray. The coincidence of these 
two signals is required for an event. This is the only appearance experiment 
and they find a signal. Once again the signal is parametrized in terms of two 
neutrino mixing. They find 

sin^ 261 - 10-3 - 3 X lO^^ (18) 
5m^ - .2 - 2 eV^ 



Summarizing neutrino oscillation data 

There are some general conclusions we may draw from the above datap-. 
The data give three distinct scales for neutrino oscillations, see figure 7Ji3 
We see 

^rnls^D » 8m\^ » Sml^i^^. (19) 

• Thus IF all three experiments are correct, we necessarily need four light 
neutrinos. This of course requires sterile neutrinos. 

Atmospheric neutrino oscillations, and perhaps even solar neutrinos, re- 
quire maximal mixing (sin'^29 ~ 1), while quark mixing angles are all small. 
We thus have 

Uc,^ » VcKM . (20) 
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Figure 6. LSND detects an anti-electron neutrino via the correlated detection of positron 
cerenkov radiation and neutron capture. 



• Hence any theory of fermion masses must give large mixing for leptons and 
small mixing for quarks. This is especially difficult in grand unified theories 
which relate quarks and leptons. 

2 "Predictive" SUSY GUT 

At this point let me define the framework for theories of fermion masses we 
call a "predictive" theory. 

1. It is a "natural" theory, i.e. if one gives you the symmetries of the theory 
and also the states and their transformations under these symmetries, 
then you can write down the Lagrangian. It is the most general La- 
grangian consistent with the symmetries. However, this is not sufficient. 

2. We also require that the number of arbitrary parameters in the theory is 
less than the number of observables. 

A predictive theory is testable. 

Clearly, to have a predictive theory one must necessarily have lots of sym- 
metry in order to reduce the number of arbitrary parameters. The appropriate 
symmetry is determined by the data. 
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Figure 7. Summary of neutrino oscillation data 



SUSY GUTs 

Consider first supersymmetric grand unification or SUSY GUTs 
SUSY dramatically restricts the form of the Lagrangian, while GUTs unify 
quarks and leptons; thus quark and lepton Yukawa matrices may be related by 
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Clebsch- Gordon coefficients. The group may be SU5, SOio or any other sim- 
ple group which breaks to the standard model gauge group at a high energy 
scale Mg- Such a theory is known to be consistent with precision clcctroweak 
data. For example, given the observed values for ai, at Mz one finds that 
the running couplings unify at a scale Mq 2 x 10^^ GeV. Running the value 
of as down from Mq to Mz one predicts the value of the strong coupling at 
the Z scale which agrees quite well with the data. 

SOio 

We use the GUT symmetry ^Oio.B This is the minimal symmetry in 
which all fermions in one family are contained in a single irreducible repre- 
sentation. We have^ 

16 D {g, u, J; /, e, v\ (21) 

In the simplest version of 5*010 one pair of Higgs doublets are contained 
in a 10. The minimal Yukawa coupling of one family of fermions to the Higgs 
fields is given by A 16 10 16 where A is the single Yukawa coupling. If we 
use this minimal coupling for the heaviest generation we obtain the symmetry 
relation Ai = Ab = Ar = A,y^ = A at the GUT scale. This is Yukawa unificatipji 
and experimentally it is known to work quite well for the third generationtS, 
but fails miserably for the other two. Using the values of mf,(m&), mT-(mT-), 
one determineSj-hpth \{Mq) and tan/3. One then predicts the value of the 
top quark masal3o 

TOt(mt) = At(mt)^ sin/3 - 170 ±20 GeV. (22) 

A brief digression: We note that the success of Yukawa unification has 
important ramifications for string theory model building. String theories, 
with SUSY GUTs realized as level one Kac-Moody algebras, necessarily need 
Wilson line breaking of the SUSY GUT down to the standard model gauge 
group. Unfortunately, Wilson line breaking preserves the nice feature of gauge 
coupling unification but typically destroys the equally nice feature of Yukawa 
coupling unification. This is a problem with this realization of GUTs in string 
theory. 

Finally, whereas £"010 unifies all quarks and leptons in a single family into 
one irreducible representation of the group, relating up to down to charged 
lepton to Dirac neutrino masses, it is not sufficient to obtain a "predictive" 

Note, V denotes a sterile neutrino. 
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theory. We need more symmetry; a family symmetry to relate states of dif- 
ferent families and to reduce the number of arbitrary parameters. Moreover 
family symmetry breaking can provide an understanding of the hierarchy of 
fermion masseSpand also ameliorate problems with flavor violating processes in 
SUSY theoriesJia To summarize, family symmetry plus GUTs are necessary 
for a "predictive" theory, but they are not necessarily sufficient. It is still 
necessary to check in each case whether or not they sufHciently constrain the 
theory. 

2.1 U2 Ui family symmetry 

In this talk, jwe consider a particular 6*010 SUSY GUT with family symmetry 
U2xUiX---M This model is a slight vatiation of the one considered previously 
in a paper by Barbieri et al. [BIIRR].EIl The changes affect the treatment of 
neutrinos without affecting the predictions for charged fermion masses. 

The three families transform as a doublet 16a, a — 1,2 and a singlet I63 
under U2 with Ui charge (-1), (0), respectively. Note, the Ui contained in U2 
counts +1 (-1) for every upper (lower) SU2 index. The hierarchy of fermion 
masses is explained by the breaking of the family symmetry. At tree level only 
the third family, I63, can couple to the 10 of Higgs bosons which also carries 
no family symmetry charge. There are three "flavon" fields in the theory 
(0'^, S'^b = S^", A°-^ = -J^") transforming (under SU2 xUiX Ui) as (2,1,0), 
(3,2,1), (1,2,2) responsible for spontaneously breaking the family symmetry. 
The vacuum expectation values [vevs] < 0^ >«< S^^ > spontaneously break 
U2 X Ui to Ui and give the second family mass. Then < A^^ > breaks the 
remaining symmetry allowing the first family to obtain massj^ 

It is important to recall that the non-abelian SU2 fanoily symmetry can 
also protect against large flavor violations such as /i — > ejS3 

2.2 The "predictive" theory 

This theory is completely defined by the symmetries, in this case 

SOwx [U2XU1X Ui{R) X Ui{PQ) X • • • ] (23) 
family symmetry 

where Ui{R), Ui{PQ) denote an R symmetry for which all fields, including 
M, M', M" have charge +1 and a Peccei-Quinn symmetry for which all 16s 
have charge +1. The ellipsis denotes additional symmetries which exist in the 

"^It is important to have only these flavon vevs in order to retain a predictive theory. 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 14 



fermion sector of the theory which may or may not remain unbroken in the 
complete theory/ 

Given the states in the theory and their charges we obtain the superspace 
potential W given by (the Ui charge is denoted by a superscript) 

— ^^charged fermion ^^neutrino ~t~ ' ' ' (^^) 

with 

Wchargedferrmon= 16° 10° 16° +16-^0° (x")' (25) 

+ {r)° [ M'° x° + 45^ 16-1] +x-^ [ M"° + 45^ 16°] 

Wneutr^no = (TE)-^^^ [Na^^^ (x")^ + N^^^ 16°] (26) 

+ (5'^^)1 Na^/^ 7V-'/^ + Na^'^ NI'^ 

The vevs of the fields Af, M\ M" give mass to the Froggatt-Nielsen0 
states labelled by x, X- The fields M', M" are SOiq singlets and we use the 
same notation for their vevs. Af, on the otherhand, is a linear superposition 
of an SOiQ singlet and adjoint with vev Mq{1 + X + (ioY) where ao, /3o 
are arbitrary constants which are fit to the data and X, Y are elements of 
the Lie algebra of SOio with X in the direction of the Ui which commutes 
with SU5 and Y the standard weak hypercharge. Each term in W has an 
arbitrary Yukawa coupling which is implicit. Finally, the fields A^o, N3 are 
SOio singlets; introduced to generate neutrino masses. These will be discussed 
in more detail shortly. 

The largest scale in the theory is assumed to be the Froggatt-Nielsen 
masses. Below this mass scale, the x, X states are integrated out of the theory 
giving the effective mass operators, see figure 8. Finally, when the "flavon" 
fields obtain vevs 

<0"> ^4>6''^ (27) 

we generate the Yukawa couplings at a scale of order the GUT scale given byQ 

e'p 
Yu = I -e'p ep reTu | A 
reTg 1 



■^One of these symmetries is needed to distinguish the fields M and M', M". 
■'^Note the parameters A, e, e', p, a, r are implicit functions of ratios of vevs, Yukawa cou- 
plings and 00 : Po- 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 15 







%-L f 10 




X X Xa f 




%-L 10 

l^M . t . t 




f Xa Xh ^ 




\ 1 Xh X 1 


1% 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 16 



Figure 8. Diagrams generating the Yukawa matrices 



Yd 
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(28) 
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1/ 



) 



A 



with 



2a 



(29) 



2cr - 1 



and 



Tf = (Baryon# - Lepton#) 
for f — {Q,u,d,L,e,v}. 



(30) 



There are six real parameters in these matrices A, e, e', p, cr, r and three 
phases which cannot be rotated away, which we take to be <I>e, <t„, $p. These 
nine parameters are then fit to the thirteen observable charged fermion masses 
and quark mixing angles. Note, once these parameters are determined the 
Dirac neutrino mass matrix rrii, is fixed. Before we consider neutrino masses 
and mixing, we should first test our theory against precision electroweak data 
which, in our case, includes the data on charged fermion masses and mixing 
angles. We shall see that the fits are quite good. 

2.3 Results for Charged Fermion Masses and Mixing Angles 

We have performed a global analysis, incorporating two (one) loop renor- 
malization group [RG] running of dimensionless (dimensionful) parameters 
from Mg to Mz in the MSSM, one loop radiative threshold corrections at 
Mz, and 3 loop QCD (1 loop QED) RG running below Mz- Electroweak 
symmetry breaking is obtained self-consistently from the effective potential 
at one loop, with all one loop threshold corrections included. This analysis 

®The predicted values of the low energy observables are highly correlated. Thus a global 
analysis is necessary in order to test the accuracy of the fit. We note that fermion masses 
and mixing angles are the precision electroweak data which constrain any theory beyond 
the Standard Model. It is important to know how well a theory beyond the Standard Model 
fits this data, even though in some cases this data still has large theoretical uncertainties. 
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Table 3. Charged fermion masses and mixing angles 

Initial parameters: 

(I/OG, Mq, ea) = (24.52, 3.05 ■ 10^^ GcV,-4.08%), 
(A,r, cr, e, p, e') = (0.79, 12.4, 0.84, 0.011, 0.043, 0.0031), 
($<T, *p) = (0.73, -1.21, 3.72)rad, 
(mo, Ml/2, M, l^{Mz)) = (1000, 300, -1437, 110) GeV, 
((m^^/mo)2, (mH„/m.o)2,tan/3) = (2.22, 1.65, 53.7) 



Observable 


Data((j) 


Theory 




(masses 


m (jeV) 


Mz 


yi.lo ( (U.uyi ) 


yi.l ( 


A/fr.r 


0\J.OOO \\J.\JO\j) 












137.04 (0.14) 


137.0 




0.1190 (0.003) 


0.1174 


Pnew ' 10 


-1.20 (1.3) 


+0.320 




173.8 (5.0) 


175.0 




4.260 (0.11) 


4.328 


Mb - Mc 


3.400 (0.2) 


3.421 




0.180 (0.050) 


0.148 


md/ms 


0.050 (0.015) 


0.0589 




0.00203 (0.00020) 


0.00201 


Mr 


1.777 (0.0018) 


1.776 




0.10566 (0.00011) 


.1057 


Me • 10^ 


0.5110 (0.00051) 


0.5110 




0.2205 (0.0026) 


0.2205 


Vcb 


0.03920 (0.0030) 


0.0403 


Vub/Vcb 


0.0800 (0.02) 


0.0691 


Bk 


0.860 (0.08) 


0.8703 


S(6^s7)-10* 


3.000 (0.47) 


2.995 


TOTAL 


3.39 





is performed using the code of Blazek et.al.. El g In this wor we present 
the results for one set of soft SUSY breaking parameters toq, M1/2 with all 
other parameters varied to obtain the best fit solution. In table || we give 
the 20 observables which enter the function, their experimental values and 

''We assume universal scalar mass mo for squarks and sleptons at Mq. We have not 
considered the flavor violating eflects of U(2) breaking scalar masses in this work. 
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the uncertainty a (in parentheses)]^ In most cases a is determined by the 1 
standard deviation experimental uncertainty, however in some cases the the- 
oretical uncertainty (~ 0.1%) inherent in our renormalization group running 
and one loop threshold corrections dominates. 

Given the 6 real Yukawa parameters and 3 complex phases we fit the 
13 fermion mass observables (charged fermion masses and mixing angles \Bk 
replacing tK as a "measure of CP violation"]); we thus have 4 predictions. 
From table ^ it is clear that this theory fits the low energy data quite well. ^ 

Finally, the squark, slepton, Higgs and gaugino spectrum of our theory 
is consistent with all available data. The lightest chargino and neutralino 
are higgsino-like with the masses close to their respective experimental limits. 
As an example of the additional predictions of this theory consider the CP 
violating mixing angles which may soon be observed at B factories. For the 
selected fit we find 

(sin 2a, sin 2/3, sin 7) = (0.74, 0.54, 0.99) (31) 

or equivalently the Wolfenstein parameters 

(p, 77) = (-0.04, 0.31) . (32) 



3 Neutrino Masses and Mixing Angles 

The parameters in the Dirac Yukawa matrix for neutrinos (eqn. ^8|) mixing 
V — V are now fixed. Of caurse, neutrino masses are much too large and we 
need to invoke the GRSY El see-saw mechanism. 



' The Jarlskog parameter J = Im(VudV*f^VcbV*j^) is a measure of CP violation. Wc test J 
by a comparison to the experimental value extracted from the well-known ii"" — A"*^ mixing 
observable ex = (2.26±0.02) X 10~^ . The largest uncertainty in such a comparison, however, 
comes in the value of the QCD bag constant Bpc- We thus exchange the Jarlskog parameter 
J for Bk in the list of low-energy data we are fitting. Our theoretical value of Bk is defined 
as that value needed to agree with ep^ for a set of fermion masses and mixing angles derived 
from the GUT-scale. We test this theoretical value against the "experimental" value of Bnq. 
This value, together with its error estimate, is obtained from recent lattice calculations.E3 
^In a future paper we intend to explore the dependence of the fits on the SUSY break- 
ing parameters and also U{2) fiavor violating effects. Note also the strange quark mass 
T7is(lGeV) ~ 150 MeV is small, consistent with recent lattice results. Finally, by adding 
one new parameter (a small ratio for the H^i to Hu Yukawa couplings) it is possible to 
obtain a small tan/3 solution that fits charged fermion masses just as well. 
'^We warn the reader that according to quite standard conventions the angle (3 is used in two 
inequivalent ways, tan/? is the ratio of Higgs vevs in the minimal supersymmetric standard 
model; while sin 2/3 refers to the CP violating angle in the unitarity triangle, measured in 
B decays. We hope that the reader can easily distinguish the two from the context. 
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Since the 16 of SO(IO) contains the "right-handed" neutrinos one pos- 
sibility is to obtain v — v Majorana masses via higher dimension operators of 
the form 

^ 16 I63 16 I63, (33) 

^ 16 163 16 i&a 

^T6 16aT6 16fc S''\ 

If 16 field gets a vev in the "right-handed" neutrino direction, then we obtain 
a, V — V mass of order < 16 >'i—JM. This possibility has been considered in 
the paper by Carone and HallB3 

The second possibility, which we follow, is to introduce SO(IO) singlet 
fields N and obtain effective mass terms v~N and N~N using only dimension 
four operators in the superspace potential. To do this, we add three new 
SO(IO) singlets {Na, a = 1,2; N3} with U(l) charges { -1/2, +1/2 }. 
These then contribute to the superspace potential (see equation Note the 
field 16 with U(l) charge —1/2 is assumed to get a vev in the "right-handed" 
neutrino direction; this vev is also needed to break the rank of SO (10). 

Finally in order to allow for the possibility of light sterile neutrinos we 
introduce a U(2) doublet of SO(IO) singlets N" or a U(2) singlet N^. They 
enter the superspace potential as follows - 

WDfi'NaN" + fisNsN^ (34) 

We show that if the dimensionful parameters /ia are of order the weak 
scale, then the sterile neutrinos (predominantly iVs) are light. The notation 
is suggestive of the similarity between these terms and the /i term in the 
Higgs sector. In both cases, we add supersymmetric mass terms and in both 
cases we need some mechanism to keep these dimensionful parameters small 
compared to the Elanck scale. This may be accomplished by symmetries, see 
for example ref. c3. 

We define the 3x3 matrix 

/// \ 

fi= \ /i' (35) 
V M3 / 

The matrix fl determines the number of coupled sterile neutrinos, i.e. there 
are 4 cases labeled by the number of neutrinos {N,y — 3, 4, 5, 6): 

• {N^ = 3) 3 active (^' = /is = 0); 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 20 



• {Ni, = 4) 3 active + 1 sterile 

ifi' = 0; fi3^ 0); 

• {Ni, = 5) 3 active + 2 sterile 

ip' ^ 0; = 0); 

• {Ni, = 6) 3 active + 3 sterile 

^0; 0); 

In this talk we consider the cases N^, — 3 and 4.0 

The generalized neutrino mass matrix is then given by F|, see figure 9 - 



where 



and 



{ 1^ N D N) 

/00m \ 

fl^ 

V 

\ ji Mn J 



7n, = Y, {Hi) = sin/3 
v2 



(36) 



(37) 



(38) 



e'Vie 

V = I -eTi6 SeVie 

re(l-a)rpl/i6 n'e, 
/O 0^ 

Mn= \qs4> 

VO0O, 

^16 J proportional to the vev of 16 (with different implicit Yukawa 

couplings) and 5", (j) are up to couplings the vevs of 5^^, 0^, respectively. 

Since both V and Mm are of order the GUT scale, the states v, N may 
be integrated out of the effective low energy theory. In this case, the effective 
neutrino mass matrix is given (at Mq) by |^ (the matrix is written in the 



This is similar to the double see-saw mechanism suggested by Mohapatra and Valle. 
™In fact, at the GUT scale Mq we define an efl'ective dimension 5 supersymmetric neutrino 
mass operator where the Higgs vev is replaced by the Higgs doublet Hu coupled to the entire 
lepton doublet. This efl'ective operator is then renormalized using one-loop renormalization 
group equations to Mz- It is only then that is replaced by its vev. 
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Figure 9. Diagrams generating the neutrino mass matrices 



{v, N) flavor basis where charged lepton masses are diagonal) 

ml^f = (39) 



'Ue 



with 

(40) 

eo = eU^ ; 1^0 = v Ul 

Ue is the 3x3 unitary matrix for left-handed leptons needed to diagonalize l^e 
(eqn. [2^ ) and cq, vq (e, v) represent the three families of left-handed leptons 
in the weak- ( mass- ) eigenstate basis for charged leptons. 

The neutrino mass matrix is diagonalized by a unitary matrix U = Uaf, 
^diag ^ jjt rnlff U* (41) 

where, in this case, a = {v,,, J^p,, i^t, i^si, v^^, Vg^} is the flavor index and 
j = {1, • • • , 6} is the neutrino mass eigenstate index. Recall, Uai is observable 
in neutrino oscillation experiments. 

In general, neutrino masses and mixing angles have many new parameters 
so that one might expect to have little predictability. However, as we shall 
now see, the U(2)xU(l) family symmetry of the theory provides a powerful 
constraint on the form of the neutrino mass matrix. In particular, the matrix 
has many zeros and few arbitrary parameters. Before discussing the four 
neutrino case, we show why 3 neutrinos cannot work without changing the 
model. 

3.1 Three neutrinos 

Consider first m^^^ for three active neutrinos. We find (at Mq) in the 
{ve, f^, Vt) basis 

'0 0' 

miff ^m' ul I & 1 I [/; (42) 
10, 



with 



, A^w^sin rri^ oj cj) 

b = ^^ + 2are 
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where in the approximation for m' we use 

V 

mti- TJitop) ^ X sin/3, (44) 

valid at the weak scale. 

rrii, is given in terms of two independent parameters { m', b }. Note, 
this theory in principle solves two problems associated with neutrino masses. 
It naturally has small mixing between — since the mixing angle comes 
purely from diagonalizing the charged lepton mass matrix which, like quarks, 
has small mixing angles. While, for b < 1, Vfj^ — v^ mixing is large without fine 
tuning. Also note, in this theory one neutrino (predominantly v,,) is massless. 

Unfortunately this theory cannot simultaneously fit both solar and atmo- 
spheric neutrino data.[^ This problem can be solved at the expense of adding 
new family symmetry breaking vevs ^ 

(5")M^(5i2)/«2-(522). (45) 

These are the most general flavor symmetry Jar eaking vevs. We discuss these 
three neutrino solutions in a future paper.Ea With ni^2 7^ the massless 
eigenvalue in the neutrino mass matrix is now lifted. This allows us to obtain a 
small mass difference between the first and second mass eigenvalues which was 
unattainable before in the large mixing limit for z/^ — i/,- . Hence a good fit to 
both solar and atmospheric neutrino data can now be found for small values of 
Ki_2. In fact, only very small values of /ci_2 are consistent with charged fermion 
masses. We note that with the addition of these two complex parameters the 
theory is no longer predictive. In the neutrino sector any of the three solar 
neutrino solutions (SMA, LMA or "Just so") can be obtained.E3 

In the next section we discuss a four neutrino solution to both solar and 
atmospheric neutrino oscillations in the theory with ki^2 = 0. 

3.2 Neutrino oscillations [ 3 active + 1 sterile ] 

In the four neutrino case the mass matrix (at Mq) is given by ^ 




(46) 



"We have checked however that we can fit both atmospheric and LSND dp 
°Such an additional vev was necessary in the analysis of Carone and Hall.| 
PThis expression defines the effective dimension 5 neutrino mass operator at Mq which is 
then renormalized to Mz in order to make contact with data. 
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where m! and h are given in eqn. ^ and 

u — a r e (47) 

V2 /13 V16 M3 Vie 

c = ~ 

LoX V sin /? Lo mt (f> 

In the analysis of neutrino masses and mixing angles we use the fits for 
charged fermion masses as input. Thus the parameter u is fixed. We then 
look for the best fit to solar and atmospheric neutrino oscillations. For this we 
use the latest Super-Kamiokande data for atmospheric neutrino oscillations u 
and the best fits to solar neutrino data including the possibility of "juat so" 
vacuum oscillations or both large and small angle MSW oscillations .□ Our 
best fit is found in tables ^ and |^. It is obtained in the following way. 

For atmospheric neutrino oscillations we have evaluated the probabili- 
ties {P{v^i v^i), P{v^ i>x) with X = {e, r, s}) as a function of 
X = Log[(L/km)/(E/GeV)]. In order to smooth out the oscillations we have 
averaged the result over a bin size, Ax — 0.5. In fig. 10a we have compared 
the results of our model with a 2 neutrino oscillation model. We see that our 
result is in good agreement with the values of 5m\^^ and sin^ ^Oatm as given. 

An approximate formula for the effective atmospheric mixing angle is 
defined by 

P{,y^ ^ i.^) = 1 - ' sin^ 20atm' sin^( " f"^ ) (48) 

4 hi 

with 

'sm^20atm' ~ 4 [ \\U^4^{1 - ||t/^4|P) (49) 

+ \\U,3f{l-\\U,s\f-\\U^4f)] 

using the approximate relation 

'Srnlt^' = Sm% « Smj^ ~ Smj^ « Smj^ ~ Sml^^. (50) 

Note, ' sin^ '29atm^ may be greater than one. This is consistent with the defini- 
tion above and also with Super-Kamiokande data where the best fit occurs for 
sin^26atm. — 1.05. In a two neutrino oscillation model, values of sin^29 > 1 
are unphysical and lead to negative probabilities. However in our four neu- 
trino model, values of ' sin^ 'iOatm'' > 1 are perfectly physical. In figure 11, we 
compare the un-averaged probability P{v^j, —f v^j) for the two neutrino model 
and the four neutrino model discussed in this talk. 

How can they give equivalent results? Experimentalists average 
the oscillating probability P{i'ij, — > v^j) over bins in the variable x = 
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♦ our model ('sin''20_ 
A (sin'2e,5mVeV') = (1.08, 3xl6'') 
▼ (sin'2e,5mVeV') = (1.08, 3.5x10'') 
X (sin 29„,^.Sm„,^VeV') = (1.08, 3.2x10'') 



■a.„,%'§<.J/eV^) = (1-08, 3.2x10"') 



1.0 
0.8 



0.6 

> 

^ 0.4 



▲ 
¥ 

T 



t 

1 



0.2 



0.0 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 

Log[(L/km)/(E/GeV)] 

Figure 10 a. Probability P(i'fi > I'/j) for atmospheric neutrinos as a function of x = 

Log[(L/km)/(E/GeV)]. This probabihty is averaged in x over bin sizes Ax = 0.5. For this 
analysis, we neglect the matter effects. 



Log[(L/km)/(E/GeV)]. This averaging process erases the unphysical neg- 
ative probabihties of the two neutrino oscillation model. It is the averaged 
result which compares well with the physical probability in a model with more 
neutrinos. Since it is much simpler to analyze the data in terms of two neu- 
trino oscillations, it is important for experimentalists to allow for unphysical 
values of sin^ 29atm in their global fits. These unphysical values may just be 
an indication of additional neutrinos. 

In fig. 10b we see however that although the atmospheric neutrino deficit 
is predominantly due to the maximal mixing between — Vt, there is nev- 
ertheless a significant 10% effect) oscillation of Vs. This effect may 
be observable at Super-Kamiokande once the ratio i?( nc jcc) (see |l^, is 
measured more accurately. 

The oscillations t'^ — > v-t or may also be visible at long baseline neutrino 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 26 



• jc = e 

♦ jc = X 
■ JC = s 



♦ ♦ ♦ 



♦ ♦ ♦ ♦ 



1.00 
0.80 
0.60 
- 0.40 
0.20 
0.00 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Log[(L/km)/(E/GeV)] 

Figure 10 b. Probabilities P{Ufj, > fx) {x = e, t and s) for atmosplieric neutrinos 



experiments. For example at K2K Ea, the mean neutrino energy E — 1.4GeV 
and distance L — 250 km corresponds to a value of x-^ 2.3 in fig. 2b and hence 
Pii^fi ^ ^'t) ^ -4 and ^ i^s) ~ .1. At Minos low energy beams with 

hybrid emulsion detectors are also being considered. These experiments can 
first test the hypothesis of muon neutrino oscillations by looking for muon 
neutrino disappearance (for x = 2.3 we have P^i^fi — > i'^) ~ .5). Verifying 
oscillations into tau or sterile neutrinos is however much more difficult. For 
example at K2K, if only quasi-elastic muon neutrino interactions (single ring 
events at SuperK) are used, then this cannot be tested. Minos, on the other 
hand, may be able to verify the oscillations into tau or sterile neutriaps by 
using the ratio of neutral current to charged current measurements (the 
so-called T test). 

For solar neutrinos we plot, in figs. 12(a,b), the probabilities {P{ve i^e), 
P{ve — > Vx) with X = {/z, r, s}) for neutrinos produced at the center of the 
sun to propagate to the surface (and then without change to earth), as a 
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Oscillation probability of v^^ to 



our model 

2 flavor model with sin '26=1.08 and 6m" = 3.2x10 eV 




■ 0.0 1.0 2.0 3.0 4.0 5.0 6.0 

Log[(L/km)/(E/GeV)] 

Figure 11. The probability P{vii — > i'^) for atmospheric neutrinos as a function of x = 
Log[(L/km)/(E/GeV)] without averaging in x. 



Table 4. Fit to atmospheric and solar neutrino 
oscillations 

Initial parameters: ( 4 neutrinos w/ large tan/3 ) 

m' = 7.11 ■ 10-2 eV , 6 = -0.521, c = 0.278, $5 = 3.40rad 



Observable 


Computed value 


SIT? 29atm 

sir? 26 sol 


3.2 • 10-^ eV^ 

1.08 
4.2 • 10-6 
3.0-10-3 



function of the neutrino energy Ej, (MeV) . |^ We compare our model to a 2 

^For this calculation we assume that electron (rie) and neutron (n„) number densities at a 
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Table 5. Neutrino Masses and Mixings 



Mass eigenvalues [eV]: 0.0, 0.002, 0.04, 0.07 
Magnitude of neutrino mixing matrix Uqi 

i = 1, • ■ ■ , 4 — labels mass eigenstates. 

a = {e, fi, T, s} labels flavor eigenstates. 



0.998 0.0204 0.0392 0.0529 
0.0689 0.291 0.567 0.767 
0.317-10-3 0.145 0.771 0.620 
0.284- 10-3 0.946 0.287 0.154 



neutrino oscillation model with the given parameters. We see that the solar 
neutrino deficit is predominantly due to the small mixing angle MSW solution 
for — Vs oscillations. The results are summarized in tables ^ and |^. 
A naive definition of the effective solar mixing angle is given by 

sin2 20i2 = 4 llt/eiiniC/eall^ (51) 

In fig. 12a we see that the naive definition of sin^ 20i2 underestimates the 
value of the effective 2 neutrino mixing angle. Thus we see that our model 
reproduces the neutrino results for Sml^i — 5m\2 = 4.2 x 10"^ eV^ but instead 
is equivalent to a 2 neutrino mixing angle sin^ 29 sol = 3 x lO^^ instead of 
sin2 26li2 = 1.6 x lO'^. 

In addition, whereas the oscillation v^ — Vs dominates we see in fig 12b that 
there is a^sigificant 8% effect) for Ve ~ v^. This result may be observable 
at SNO with threshold E>h MeV for which P(ye Vf,) ~ .05. 

These results may^be compared with a recent analysis by Bahcall, Krastev 
and Smirnov [BKS98P, see figure 13. The shaded area is the 99% CL fit 
for sterile neutrinos consistent with the total rates, zenith-angle and recoil 
electron energy spectrum. The large dot is the best fit and the cross is our 
solution. 

We note that, even though we have four neutrinos, we are not able to 
simultaneously fit atmospheric, solar and LSND data, i.e. it is not possible to 
get both "(5m^^_j^^" and sin^20 large enough to be consistent with LSND. It is 
worth asking however how robust is this result. Not surprisingly we find that 

distance r from the center of the sun are given by (rie, n„) = (4.6, 2.2) X 10^^ exp(— 10.5-^) 
eV^ where i? is a solar radius. We also use an analytic approximation necessary to accpunt 
for both large and small oscillation scales. For the details, see the forthcoming paper .^j 
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♦ ♦ our model (sin'lGjjjSmVeV^) = (1.6xl0"\ 4.2x10") 

▼ ▼ (sin'2e,5mVeV') = (2x10 , 4.2x10'') 

A ▲ (sin'2e,5mVeV') = (8x10"', 4.2x10'') 

• • (sin'29,„,.Sm,„,'/eV') = (3x10"', 4.2x10'') 




(MeV) 

Figure 12 a. Probability P{ue > Ue) for solar neutrinos 



upon introducing the new parameters K1.2 (eqn. ^5|) a^olution to all three - 
atmospheric, solar and LSND oscillations is obtained.E^ 

Finally let's discuss whether the parameters necessary for the fit make 
sense. We have three arbitrary parameters. We have taken b complex, while 
any phases for m' and c are unobservable. A large mixing angle for i/^ — v^- 
oscillations is obtained with |6| ^ 0.5. This does not require any fine tuning; 
it is consistent with -^-^^ ~ 1 which is perfectly natural (see eqn. ^3|). The 

parameter c[eqn.|^ « 0.28 « implies ^13 w 0.41 ^ mf. Thus in 

order to have a light sterile neutrino we need the parameter fi^ ~ 70 GeV 
for (p ~ Vig. Considering that the standard fi parameter (see the parameter 
list in the captions to table |^) with value /i = 110 GeV and /is [eqn. Q 
may have similar origins, both generated once SUSY is spontaneously broken, 
we feel that it is natural to have a light sterile neutrino. Lastly consider the 
overall scale of symmetry breaking, i.e. the see-saw scale. We have m' = 



ws-p8-50x6-00: submitted to World Scientific on February 1, 2008 30 




(MeV) 

Figure 12 b. Probabilities P{i'e > Ux) (x = /i, t and s) for solar neutrinos 



7 X 10-2eV[table | « Thus we find ^^^^ - ^ - 6 x IQi* GeV. 

This is perfectly reasonable for (16) ^ (0^) ~ Mq once the implicit Yukawa 
couplings are taken into account. 

4 Conclusions and future tests 

We have discussed "predictive" theories of charged fermion masses in this talk. 
These theories are constrained by lots of symmetry, in particular SUSY grand 
unification and family symmetries. We have argued that such "predictive" 
theories may constrain neutrino masses and mixing angles as well. And vice 
versa, neutrino data will therefore be able to constrain theories of charged 
fermion masses. 

We discussed a particular model with an SOiq x C/2 x C^i x • ■ ■ symmetry. 
With minimal family symmetry breaking vevs of the form (fp ^ S^"^ > A^^ this 
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4T 



Figure 13. Fit to solar neutrino oscillations by BKS 98. 



theory is predictive. It leads to a very simple picture for neutrino oscillations. 
We find maximal j/^ — s- Vr mixing for atmospheric neutrino oscillations; a 
small mixing angle MSW solution with ^ i^steriie mixing for solar oscilla- 
tions and NO solution for LSND. 

We also made the important point that the origin of light sterile neutrinos 
in SUSY theories may be related to the origin of the Higgs /i term. In order 
to obtain one light sterile neutrino in our model we needed a supersymmetric 
mass term of the form ^^N^N^ where both N3, are SOio singlets and 113 
is of order /i. 

Our model will be tested in future neutrino experiments. 

• For atmospheric neutrino oscillations we have — s- Super- 
Kamiokande will be able to distinguish between and Vsteriie by look- 
ing at the ratio of neutral current to charged current processes. A re- 
cent analysis by Super-KamiokandeQ measures the ratio R[nc/cc) = 
, Z'^)-""'" — consistent with one with large errors. This result favors 
Vr but the result is not yet significant. In addition, by looking at the 
zenith-angle dependence one may be able to distinguish between Vr , i^s 
as well sincCpfpr Vs, but npt I'r, there is an MSW effect in the earth. 
Finally, K2KE3 and MINOSHj wiU both be able to confirm the disappear- 
ance of j/^ and MINOS should eventually be able to see the appearance 

of Vt- 

• For solar neutrino oscillations we find a snmJJ mixing angle MSW solution 
with ^ ^s- Both Super-K and BorexinoEj should be able to distinguish 
between small mixing angle, large mixing angle or vacuum oscillation set: 
lutions by their characteristic energy dependence and seasonal variation^ 
In addition, SNO should be able to distinguish between active and sterile 
neutrinos by measuring the ratio of neutral current to charged current 
solar neutrino processes. A ratio of order one is indicative of i^e ~* inactive, 
while a ratio much less than one would confirm Ve — > Vsteriie oscillations. 
Finally, KAMLANDEI (similar to the CHOOZ experiment) will be sen- 
sitive to the large angle MSW oscillation region. It can thus confirm or 
rule out this possibility. 

• We find no evidence foriSND oscillations. Fortunately, this result will be 
tested by MiniB00NEc3; a short base-line oscillation experiment which 
will be done at Fermilab. 
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